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OS ' Abstract In [FMRT1] we introduced a multiscale expansion for many Fermion systems in two 

space dimensions based on a so-called sector decomposition. In this paper a completely different 
expansion is introduced to treat the more difficult case of three (or more) space dimensions: it is 
' based on an auxiliary scale decomposition and the use of the Hadamard inequality. We prove that 

^| . the perturbative expansion for a single scale model has a convergence radius independent of the 

scale. This is a typical result, already proved in the two dimensional case in [FMRT1], which we 
cannot obtain in three dimensions by naive extrapolation of the sector method. Although we do 
not treat in this paper the full (multiscale) system, we hope this new method to be a significant 
step towards the rigorous construction of the BCS theory of superconductivity in three space 
i dimensions. 
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I Introduction 



In this paper we consider many Fermion systems formally characterized by the effective po- 



£H ■ tential 

g ; = io g \\ e - v ^ + ^+^ c w>,vo (i.o) 



for the external fields ip e , i\) e . Here, d\xc{^, VO is the Fermionic Gaussian measure in the Grassmann 
variables {tp(£), ?/>(£) | £ G IR x IR rf x {f , J.}} with propagator 

f d d+1 n p i <P>?^?>- 

c &£) = S °>* / TT ( L1 ) 

J (27r) d+i ip - e(p) 

where 

<p,(, >-=p-x-p T (1.2) 

and 

e(p) = |--M- (1-3) 

The variable £ = (r, x, a) consists of time, space and spin components and the (d + l)-momentum 
p = (po, p). The interaction is given by 

A r 4 

V= ^jU^ ^^2,6^4)^(^)^(6)^(^)^(6) (1-4) 



where A is the coupling constant, the kernel V(£i, £2, £3, £4) is translation invariant with V(0, £2, £3, £4) 
integrable and 

J ~T~,-.JtR JTR d 



*€{U}' 



The partition function is 



Z = J e- y ^^cW>,V0 (1.6) 
so that 0(0,0) = 0. 

The fact that V is non-local (although short range) is the source of a few inessential complica- 
tions. Therefore in this paper we restrict ourselves to the case of a 5 function to focus the reader's 
attention on essential aspects. Therefore from now on, we will consider 

Va = ^ / dxdT ( x > T ) i>{ ( x > T ) ( x > T )^l ( x > r ) ( IA - b ) 



2 J A 

The Euclidean Green's functions with 2p points 

p S 2 

Opto, £1, q = n ^ e(6)(5lM / (i-7) 

generated by the effective potential are the connected Green's functions amputated by the free 
propagator. By definition, Q exists when the norm 

||G P || =max sup / TTd& |G p (£i, • • • , & P )| (1-8) 

of each of its moments, G p , p > 1, is finite. Intuitively, ||G P || is the supremum in momentum space 
of G p . In fact, the supremum in momentum space was used as the standard norm on vertices in 
[FT2] . 

Our goal is to give a rigorous proof that the standard model for a weakly interacting system 
of electrons and phonons has a superconducting ground state at sufficiently low temperature. 
Perturbation theory and, in particular, the renormalization of the two point function was controlled 
in [FT1]. A renormalization group flow for the four point function was defined and analyzed in 
[FT2]. Two additional ingredients are required to complete this program. First we need an infinite 
volume expansion that controls non-perturbatively this renormalization group flow. This first 
ingredient, "constructive stability" was provided in two space dimensions (d=2) in [FMRT1]. The 
other physically interesting case, namely d = 3, is much harder. In this paper we prove by a 
completely different method that one of the key results of [FMRT1] extends to d = 3: the radius 
of convergence of the single scale theory is uniformly bounded below. However up to now, we have 
not been able to extend this result to the multiscale theory. Therefore "constructive stability" in 
d = 3 remains an open problem. 

The second ingredient is the control of the reduced BCS model for Cooper pairs via another 
expansion (of the 1/N type). One has to prove existence of the BCS gap, of the associated 
spontaneous 17(1) symmetry breaking, and of the masslessness of the associated Goldstone boson. 
This program is well under way. For references on this part of the program and on the general 
strategy see [FMRT3-5] . 

In the remainder of this paper we specialize to d = 3. Our method extends however without 
difficulty to any d > 2 (see Appendix I). 



2 



As in [FT1,2], the model is sliced into energy regimes by decomposing momentum space into 
shells around the Fermi surface. The j th slice has covariancc 

where 

/ J -(p) = /(M- ! «(pg + e(p) 2 )) (1.10) 

effectively forces \ipo — e(p)| ~ MK The function / £ Cq°([1,M 4 ]). The parameter M is strictly 
bigger than one so that the scales near the Fermi surface have j near — oo. The model is defined 
in finite volume and at fixed scale by the following lemma: 



Lemma 1.1 



where 



A f 4 

VA =2y A4 n^ ^i,&,&,€4) ^1)^(6)^4)^3) (1.12) 



and 



7 U) 



J e- v M'®dn cU) (il>,ii) (1.13) 

is analytic in X in a neighborhood of the origin that includes the disk of radius const (M 2j '|A|) 1 . 

The full proof (not very difficult) is given in [FMRT1, Lemma 1]. Remark however that the 
radius of convergence depends on volume and scale in a patently unsatisfactory way. 

Recall that Gp' A ^ is the 2p-point Green's function generated by Q^} ■ By the Lemma above 
the Taylor series 

00 

=J29n(p,j,A)X n (L14) 

ra=0 

has a strictly positive, though possibly j and A dependent, radius of convergence. The main result 
of this paper is 

Theorem I There exists a const, independent of j and A, such that 

\\g n (p,j,A)\\ < const "+PM< 4 - 2 ^ II F|r (1.15) 

where 

\\9n(p,j, A) || = max s 
Furthermore the limits 

9n(p,j)= lim g n (p,j,A) (1.17) 

A^IR 3 

exist and the infinite volume Green's functions at scale j 

00 

f-EfctPjlA" (1.18) 

n=0 

are analytic in |A| < R = (const ||V||) _1 . 



g n (p,j,A)\\ =max sup / TT d& \g n (p, j, A)(£i, • • • , £ 2p )|. (1.16) 



This theorem is analogous to the first theorem of [FMRT1], which applied only to d = 2. Only 
the dependence in the external fields is different. Remark that this external field dependence, i.e. 
the factor M^ 4 " 2 ^ is still not the desired perturbative one which would be M 1 / 2 ^ 2 ^' [FT1]. 

A theorem such as Theorem I is usually proven with a standard cluster/Mayer expansion. 
Space-time, lR d+1 , is paved by cubes A of side M -3 dual to the decay rate M J of the propagator. 
The decay rate is primarily determined by the thickness of the shell in momentum space. Then 
one expands in coupling constants that control the interaction between boxes. One essential 
prerequisite for the convergence of such expansions is the j independent estimate |Z(A)| < const . 
For our models, one can see in perturbation theory that this estimate fails. For instance the 
vacuum graph G2 (see Fig. I) can be computed explicitly and in d = 3 it leads to a second order 
contribution in e^^ 3 to \Z(A)\. 

This difficulty can be traced to the fact that the Pauli exclusion principle now permits about 
M - ( d-1 W electrons to be located in A with momentum restricted to the shell. For d = 1 the Fermi 
surface consists of just two points and consequently only O(l) electrons are allowed in A and there 
is no difficulty. As d grows the Pauli exclusion principle becomes progressively weaker and the 
estimate on the partition function in A becomes more and more j dependent. 

In [FMRT1] the solution to this difficulty was found in a decomposition of the shell into smaller 
units, called sectors. This idea does not seem to work beyond two dimensions. An informal 
discussion of why this is so is included for the interested reader in Appendix II. Our proof of 
Theorem I relies therefore on a completely different idea : we decompose the propagator in a slice 
according to an auxiliary scale, and the Hadamard inequality is used to control the Fermionic 
determinants. 

Acknowledgements This paper is part of a larger program in collaboration with J. Feldman 
and E. Trubowitz. We thank particularly J. Feldman for his comments. 

II The auxiliary scale decomposition 



The j-th slice propagator is 



J (2tt) 4 ip - e(p) 

1 f +oc , lP0 (f-T) [ +OC , > n , .47rsimymi|x-x| 
— j / d Po e'P^ T / du f j (p ,u)D( Po ,u) K — J: l - 

(27TJ 4 J-00 J-i 2|x-x| 

+M j+2 r +M 3+2 



= 5**-^, -/ / dpodu e ipo{f - T) f j (p ,u)D(p Q ,u)sm(VT+^\ X - (III) 

by putting u = e(p), fj(po,u) = f{M~ 2 3{p\ +u 2 )), and D(po,u) = (ipo — u)^ 1 . This formula is 
proved by performing the integrals over the angular components 9 and of p. In the last line we 
used the fact that / has its support in [1, M 4 ], and we assumed that j < —2*. 

We decompose further this propagator into auxiliary slices according to the scale of |x — x| by 
writing 



cHt,0 = Y,c J ' k (^0 (n.2) 

k=j 



* The first slices j = or j = — 1 are trivial from the point of view of this paper. 
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with: 

87HX - X| 7-M3+2 J-M3 + 2 

D(p ,u)fj(p , u) sin(Vl + u|x - x|) (II. 3) 

^°(x,x) = e-l x ~*l 2 (II.4a) 
^ fc (x,x) = e - M2fc l x -"l 2 - e -M 2fc+2 |x-x| 2 for - 1 > k > j + 1 (II.4b) 

A** (x, X) = 1 - e -M 2j + 2 |x-x| 2 (H.4c) 

These formulas are totally explicit. We choose them so that C- 7 ' is a nice convolution of C- 7 with 
an explicit spatial Gaussian, since the Fourier transform of a Gaussian is Gaussian, but other 
formulas are possible. We have 

Lemma II. 1 There exists a constant K(p) such that, for k = 0,...,j, and any integer p the 
functions C^ k obey the bounds : 

|C J ' fe (C,«e)| < K{p)M j+k {l + M k \x-x\)- p {l + M j \r - f\)~ p (II.5a) 

Furthermore 

C J '' fe (£,0 = for —l>k>j; |C i,0 (£,0| < M 2j (II.5b) 



Proof : For k ^ j, we use the identity 



(1 + iM 3 <r - f)fe lpa{f - T) = (1 + M j 4-) P e iMf ~ T) 

dpo 

and the obvious inequality (1 + M J |r — f |) < \/2\l + iM^(t — f )| to write 

r+M j + 2 r +M j+2 



(11.6) 



(1 + AP\t -f\) p \C^ k \ < (V2) p 



^ /—\-lVl" /-fM" 

( (1 + Mi ^ PeiMf " T) ) ^ (:P °' U)/ ^ P0 ' U) sin (\/l+^|x - x| 



dp Q du 



= (V2) P 



^< fc (x,x)<W — / / dp due ip °^ 

87T d |x - X| J_ M 3+2 J-M1+ 2 

( (1 " Mj ^) PD (P°' u )fj(Po> u )) sin(>/TTS|x-x|) 



(II.7) 



The last line is obtained by integration by parts; boundary terms vanish by the support properties 
of fj. Computing the action of partial derivatives M- 7 ^ on D(po,u)fj(p ,u) completes the proof. 
Indeed by the scaling and support properties of fj we obtain, up to a constant K(p), the same bound 
as for the initial function D(p , u)fj(p , u). More precisely for k ^ we bound sin(\/l + u\x — x|) 
by 1. Using then 

f -+M J + 2 r +M j+2 

(Hi 



we obtain that 



/ + iW J 1 r-\-M* 
/ dp du\D(p ,u)\ < M j . 

-Mi+ 2 J-M3+ 2 

Ai' k (x,x) 



(l + M J |r-f|) p |C J ' fc | < ^i(p)M J 



X — X 



(II.9) 
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and we remark that for k ^ j, k ^ and for any p 



x — x 



< K 2 {p)M k {\ + M fc |x - x|)- p 



(11.10) 



(In fact the function A k has exponential decrease for k ^ j). In the case k = 0, we use the bound 

sin(v^|x-x|) ^ 
|x — x| 

for |m| < 1 (which is true on the integration domain). We combine this bound with the trivial 
estimate 



A^°(x,x) 



<*r 2 (p)(i + |x-x|)- 



(11.12) 



to obtain (II. 5a). 

In the last case k = j, the spatial decay also has to come from integration by parts. For 
simplicity we assume that p = 2q is even (this obviously implies the general case as well) . We use 
the identity: 

M 2j '|x-x| 2 sin(VTT^|x-x|) = (2iM 3 VTTTi-^) 2 sin(\/TT^|x - x|) (11.13) 

du 

Therefore we have 



(1 + M J '|x - x|) p (l + M 3 \t - f|) p |C"'| < 2 P 



(1 + M 2j |x - x| 2 )«(l + iM\t - f)) p C J " 



2" 



-/ / dp du[{l + M ] —fe tMf - T) 



87T 3 |X 



dp ' 



D(p Q , u)fj(po, u) ((1 + (2zM^VTT^) 2 ) 9 sin(yTT^|x - x|)^ 



(11.14) 



Then we integrate by parts, both on p and u. The boundary terms vanish again by the support 
properties of fj, hence partial derivatives (1 — M 3 -^) act as before on D(po,u)fj(po,u), and 
operators 2iM 3 \J\ + u ^ change into transposed operators (— 2iM 3 -^£)\/\ + u ■ , also acting on 
D(po,u)fj(po,u). This action, by Leibniz rule, is slightly complicated to write down. But clearly, 
again by the support and scaling properties of fj, we get up to a constant K(p) the same estimate 
than for the initial po and u integral, namely M 3 . We bound | sin(\/l + u|x — x|)| by one, and the 
last factor M k = M 3 for k = j comes from the easy estimate on (II. 4c): 



A^(x,x) 



< M 3 



(11.15) 



Finally (II. 5b) is easy and left to the reader. (For the case k = 0, use the fact that 
/ dp duD(po,u)fj(p ,u) = 0.) ■ 
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Corresponding to this division of the covariance, there is an associated orthogonal decom- 
position of the fields ^ and $ distributed with covariance as sums of fields and 
distributed with covariance C J ' fe : 



For each scale k we shall perform in the next section a cluster expansion according to a lattice 
T>k of tubes of spatial side M~ h and of time side M - - 7 . The power counting with respect to this 
auxiliary scale is superrenormalizable and identical to the one of the ultraviolet limit of the <^|. 
To understand this fact, we can compute the weight of a connected graph with n vertices and I 
internal lines at subscale (j, k) with k > j and e = 2p = An — 21 external legs. Using the same norm 
as in Theorem I of section 1, we fix a particular vertex and integrate over the others. The usual 
estimate, selecting a tree in the graph [R], gives M _ ( J ' +3fc ^ n_1 ) for integration of the vertices times 
M^+W for the power counting of the lines. This is M^-^W^-^Xi"''). The factor M^ 4 " 2 ^ 
is that of (1.15) in Theorem I. The sum over k is convergent provided n + p — 3 > 0, which selects 
exactly the same graphs as the condition of ultraviolet convergence in the <j>\ theory. This is to 
our knowledge the first example of a Fermionic model with this kind of power counting. 



We conclude that as far as renormalization is concerned, the single slice model treated in this 
paper is similar in difficulty to the old ultraviolet construction of by now well established [G J] . 
Vacuum graphs are divergent up to three vertices hence the graphs G\ G 2 , G' 2 and G3, G' 3 , G' 3 ... 
of Fig. 1, are divergent and the two point subgraphs are divergent up to two vertices, hence the 
graphs B\ and B 2 , B' 2 of Fig. 1, are a priori divergent. But by (II. 5b) the "tadpole" is convergent, 
hence G\, B\ and all the "primed" graphs in fact converge. The second order vacuum graph G 2 
is linearly divergent and has to be cancelled explicitly by going to intensive quantities such as the 
pressure or normalized Schwinger functions. Finally the graph B 2 is apparently logarithmically 
divergent if we replace the propagators by the bound (II. 5a). But in fact it is convergent because 
of the oscillations of the sin 3 function, which are lost in the right hand side of (II. 5a). To exploit 
this fact, one must renormalize this graph, hence make an explicit computation for it. Had this 
graph not been convergent, we would have been able only to prove a bound in |j| -1 / 2 on the radius 
of convergence in Theorem I. 

Ill The expansion 

Since this problem is superrenormalizable, we can treat the few divergent graphs by introducing 
explicit counterterms, in the manner of [G J] . This is no problem for the vacuum graphs since they 
will quotient out anyway in the normalized functions. But for the two point subgraphs B\ and B 2 , 
we have to compensate for these counterterms by changing the propagator: 

Lemma III.l The infinite volume connected Green's functions defined by 



(0 



(11.16) 



Figure 1: The divergent graphs in d=3 



p 



Gp(£i,£i 



Q 



(III.1) 



1=1 




(HI.2) 



are identical to those of the theory 




(III.3) 
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where 

oo 

C(p) = J2C(p)(B(\,p)C(p)) n (III.4) 

n=0 

1 = V(tp + ^ e ,1p + tpe) + M + W (111.5) 

M = xJ ^5^(0^(0 + A 2 J d^db $(£i)B 2 (£i, (IH.6) 

W = J (a\ + b\ 2 + c\ 3 ) (111.7) 
where B(X,p) is the Fourier transform of the translation invariant kernel XBi5(£i, &)+X 2 B 2 (t;i, £2)- 



Proof : In the perturbation expansion the counterterm W cancels out, and the formal power 
series expansion of the propagator C exactly cancels against the insertions of the mass counterterm 
M to give back the original propagator C. But by Lemma 1.1, the perturbation expansion really 
defines the theory. This proves that the two theories are also identical at the constructive level, as 
defined by the same analytic germ. ■ 

We apply this lemma with a b and c in W the amplitudes of the divergent vacuum graphs 
with one, two and three vertices in Fig. 1, and 

M = xj ^1^(0^(0 + a 2 / dkdfc <Ka)52(a,6M&) ■ (in.8) 

where B\ and B 2 (p) are the amplitudes for the corresponding graphs in Fig. 1. The tadpole 
amplitude B\ is 

Si = ci / dq - /j(<?) , , (HL9) 
J iq - e(q) 

The Fourier transform B 2 (p) of the kernel B2 is exactly 

nt\ [j a fj(P + Qi) fj(-Qi+Q2) fj(Q2) mnn > 

B 2 (p) = c 2 / d qi dq 2 -— J - — ■ ■ — — ^ — - (IIL10 

J i(P + qi)o - e(p + qi) i(-qi + q 2 )o - e(-qi + q 2 ) i(q 2 )o - e(q 2 ) 

(ci and c 2 are some inessential numerical constant). 

Lemma III. 2 For p in the support of ' C J , there exists some numerical constant K such that 

\B X \<KM^ (III.ll) 
\B 2 (p)\ < KM 2j (111.12) 



Proof : (III.ll) follows from (II. 5b). To bound B 2 we can bound the three denominators in 
(III. 10) by M~ 3j , the two integrals on dqifl and dq 2 $ by M 2j , and the two spatial integrals on 
dqi and dq 2 by M 3j . Indeed let S a be the sphere of radius 1 centered around a. The qi integral 
is restricted to distance M- 7 of S- p by the function e(p + qi), hence runs on a volume of size M 3 . 
But for fixed qi generic, i.e. such that |qi| is not close to or 2 v / 2m/x, the q 2 integral is restricted 
to distance ~ M 3 of the circle So H 5_ qi , hence runs on a volume of size 0(l)M 2j . The degenerate 
cases |qi| — 2yJ2mfj, leads also to the same bound. ■ 
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Remark also that B 2 {p) is spatially rotation invariant so that it can be written solely as a 
function of po and u: 

B 2 {p) = B 2 (p ,u) (III. 13) 

with u = e(p). 

By lemma III. 2, the power series &{p) — X^^Lo ^ j '(p) ^(•^i + ^ 2 B 2 (p))C :l (p)^j are conver- 
gent, uniformly in j and p, for A < 0(1)*. 

The propagator & is given by an integral representation identical to (II. 1) but with D(p , u) = 
(ipo — u)^ 1 replaced by D(p ,u) = (ipo — u — fj{po,u)(\B\ + A 2 B 2 (po, «))) We can slice the 
propagator C' J = *Yl,kC^ k exactly as in (II. 3). Lemma 1 then applies exactly in the same way, 
since D and its derivatives satisfy the same bounds as D on the support of fj. This proves: 

Lemma III. 3 The sliced propagator &' satisfies the same bounds (II.5a-b) as C^ k . 

From now on we will forget the bars, and for simplicity we write C- 7 , C J ' fe , D instead of & 
&' k , D. 

To obtain a better radius of convergence for our single slice model, we return to the finite 
volume theory. For this theory, the interaction 2a is equal to T restricted to A, hence is: 



Ia = V a + M x + W a 

a r 4 

Va= 2/^11^ ^(£i>6>6>£0^(W(W(W(6) (ni.14) 

2 

m a = a / dtBMom + [ ^(a)£ 2 (a, (m.15) 
j a j a- t\ 

Wa = |A|(aA + 6A 2 + cA 3 ) (111.16) 

where Wa stands for the finite volume vacuum counterterms, and Ai\ for the finite volume two 
point function counterterms. 

Remark that the finite volume Green's functions for that theory this time differ by a (small) 
change in the boundary condition from those of the initial theory put in a finite volume. This 
change in the boundary condition comes from the fact that the counterterms are limited to the 
volume A, whether the power series for the propagator is computed as if they where introduced in 
the whole space H d+1 . But by lemma 1.1 we know that we are in a single phase, at least for A small 
enough, and that the thermodynamic limit is independent of boundary conditions. Therefore we 
can compute it in this way as well. 

The theory is analyzed as usual with a sequence of cluster expansions. Let T>j be a covering 
of TR d+1 by hypercubes of side size M _J , and 2?fc be a covering of IR d+1 by hyperrectangles which 
have size M~ k in each spatial direction and size M^^ in the (imaginary) time direction. 

The finite volume A is chosen to be a finite union of AT (A) hypercubes in T>j. This union is 
called Vj(A). We choose M integer so that each such hypercube is then decomposed into M -3 ( J_fe ) 
hyperrectangles of T>k ■ Therefore A is also a union T>k (A) of hyperrectangles of T>k ■ 

* By a fixed point analysis we could iterate this argument so that B\ and B 2 themselves 
are replaced by similar quantities called _Biand B 2 , computed with propagators C J instead of 
C J '(this is the so-called problem of "bubbles within bubbles"). This is not really necessary for our 
construction, but would simplify the cancellation of the graphs B\ and B 2 with the counterterms 
of M. 
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For A € T>k we call \A the characteristic function of A. We decompose the external functional 
derivatives in Q according to the cubes of T>j(A) in which they act (recall that one of them is fixed, 
say at the origin, and the others have to be summed in A) . For any subset X of hyperrectanglcs 
at higher scales we can also define Qx as those derivatives but restricted to the support of X. 
We perform the first horizontal cluster expansion with respect to the k — scale. Let us denote 
by <><fe, for k = 0, the (normalized) Grassmann integral with respect to the Gaussian 
(anticommuting) measure of propagator XZ;<fc • 

We consider a given Green's function 

(G P )< = -±-J Gpe-^^d^ «>, VO , (IH-17) 
where Z A = J e-^^^dncu) W>, i>) and 

2p 2p 



G p 



U^j G(£i,---,6 P )I1$(&) ( IIL18 ) 

J=l 3=1 



is an arbitrary monomial. 

The non-local countertcrm B2 is decomposed as 



B 2 {s,o = c(y2c j *{t,o) (ni.19) 
^ k ' 

A standard cluster expansion with respect to C- 7,0 is now performed, including the propagators in 
the non-local counterterm B 2 - The formula is ([B][AR]) : 

(G)< = ( E Wifl A w) (ni.20) 



Q 

H i=l 



The polymers Xi are non-empty sets of hyperrectangles of 2?0i an d the amplitude A(Xi) for the 
polymer Xi is given by a "Brydges formula" 

= Y ^J dTh J ^^(^{n^jfoe-^^* (IH-21) 

where Zq is the normalization factor of a trivial polymer made of a single hyperrectangle with no 
external legs. The sum is performed over oriented trees built on the polymer X. The orientation 
tells us when a tree line I joins X to X' if it is £ which is located in X and £ in X', in which case 
we put A/ = A and A; = A', or the contrary, in which case we put A; = A and A; = A'. The 
linking operators are 

Li= [ £ ci,0 (z>z) xl s ( ~ SJ s ,* 

J A, Ja 1 <tyo(0 <ty>0(£) 

+3 C / d£ / d^(^(o^''°(^of^ o (T,/ l ,e,o+ E ^' fc (^o) ( IIL22 ) 

Jai Ja, \ k<-\ ' 

The second term corresponds to the non-local coupling through the B2 counterterm. 

d T h is the ordinary Lebesgue measure over interpolation parameters h = {hi,h <G T}, each 
running in [0,1]. dfi C i-°(T.h) 1S the normalized Grassmann Gaussian measure with propagator 
C j ' (T,h), depending on the interpolation parameters h, which is equal to C J '°(^,^) ■ ^ T (£,£), 
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where ft- T (£,£) is if £ or £ are not in the support X of the polymer, is 1 if they lie in the same 
hyperrectangle of X, and is otherwise the infimum over all parameters hi on the unique path 
joining the hyperrectangle containing £ to the hyperrectangle containing £ [B][AR]. 

Then we apply a vertical expansion which tests whether a polymer X contains or not fields of 
lower scales (recall that external fields in Q are considered by convention of the scale j — 1 , so lower 
than all other fields). This means that we introduce a parameter v x which for each X interpolates: 

Vx = V x + + V^ linkod to (111.23) 

V x = \f dxdT^?(x,r)^(x,r)^(x,r)^(x,r) (III.24a) 

Vf" 1 = \ ^^- 1 (x,r)^f- 1 (x,r)^- 1 (x,T)Vf- 1 (x,r) (III.24b) 

y linked to <-l = 1 ^ rfx ^ T(X)T) ^ (XjT)V , T(X)T)V , i(X)r) 

(III.24c) 

Vx(vx) =V° X + V x - 1 + v x V° x linkcd to ^ ; V X = V x (v x )\ vx=1 (111.25) 

and we also interpolate M. x in a similar way. Then we perform a Taylor expansion in v x to first 
order for each X. The terms with v x = are called vacuum polymers. For the error terms we 
draw a link between the polymer X and all the hyperrectangles of the next scale containing a 
hyperrectangle of its support; in this way we construct a polymer X living on the hyperrectangles 
of the next scale. 

Then we apply a contraction rule (integration by parts on the field) , to compensate explicitly 
the vacuum polymers which contain up to three vertices with the corresponding counterterms, and 
the polymers with exactly one or two vertices and two low momentum fields of scale < — 1. This 
is standard [GJ]. There remains only polymers containing at least 4 vertices or at least 4 external 
legs (a B 2 counterterm counts for 2 vertices, etc.). 

This explicit rule a la Glimm-Jaffe is slightly simpler than the modern more powerful renor- 
malization group schemes that are needed for just renormalizable theories [B],[R]. Indeed it avoids 
to perform at each scale a Mayer expansion to factorize all the vacuum graphs, and e.g. exponen- 
tiate all the two point functions. The result of a Mayer expansion is a formula expressed in terms 
of "Mayer configurations" , i.e. finite sequences of polymers. Iteration of Mayer expansions lead to 
sequences of sequences, hence to formulas heavier to manipulate than ordinary polymers. 

To iterate we apply the second scale cluster expansion to the set of non trivial subsets of D_i 
whose elements are the polymers obtained after the first expansion, and we apply a second cluster 
expansion with respect to the measure C 3 ' -1 between these units, etc... 

We obtain the following final formula 

(s)< = E hfi A w ( IIL27 ) 

Yl,-,Yq q ' i=l 

Yf(Y),- -,Y£(Y) disjoint 

The polymers Y at the last scale have vertical tree structure, that is they are a collection of 1 + \j \ 
sets Y k for k = 0, • • • , j. More precisely, Y k (Y) is a subset of T>k, i.e. a set of hyperrectangles, 
made of connected components Y t k (taking at level k into account the connections of scales higher 
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than or equal to k, and the set of all the has tree structure for the inclusion relation (see e.g. 
[R]). 

The amplitude A(Yi) for the polymer Yi is given by applying inductively the Brydges formula, 
according to the collection of subsets X k (Y). 

A(Y)= ]T ...fd T s f d^J^L'Xgxe-^^ (111.28) 

T ,...,Tj J J ^l£T ' 

where the links operators L\ include the explicit contraction rules to cancel the counterterms for 
the graphs B\, B 2 and Gi,..., G3 and the sum in (III. 28) is restricted by the condition that any 
connected component at a given level has either more than three vertices or more than two external 
legs. 

Theorem I then follows from the following bound 
Theorem III.l For any K > there exists some r{K) > such that for |A| < r(K) 

\A(Y)\K^ X ° (Y ^ < 1 (111.29) 

Y\0£X°(Y) 

Indeed performing a standard Mayer expansion we obtain analyticity of < Q > for |A| < r(e) 
(where e = 2.718..). 

The next section is devoted to a proof of Theorem III.l. 

IV The bounds 

The formula giving the amplitude of a polymer Y is expanded into a sum of diagrams made 
of explicit propagators derived by the cluster expansion, vertices, and fields. The Fermionic Gaus- 
sian integration on the fields is simply a formal device which changes the fields into products of 
determinants. We recall the "method of combinatorial factors" to keep track of the many sums in 
the expansion. This technique uses the elementary estimate 



K, > , 2J K 4 1 < 1 



<sap\KiUi\. (IV.l) 



to replace each sum by a suprcmum. To help remember the combinatorial factor Ki multiplying 
the value U{ of a given diagram, the factor will be assigned to a specific line or vertex or field of 
the diagram. Remember also that in the end we have a small coupling constant per vertex, hence 
per field, so that such constants can be forgotten in the estimates. 

A vertex will be localized in the tube of its highest leg. That means that we write 

r H fc < A > _ H fc l <-) fc 2 (-) fc 3 

e v Y=e E &€Y f A ^ (0E» 1 . fca ,*,< M A,^ («) 
II E (/ «^ A) «5 £ ^^ k2 (0 ( f\0) (iv.2) 

A „(A)=0 1 h V A fei,fe 2 ,fe3<fe(A) / 

Furthermore by the usual "local factorial principle" we can extract, up to a constant per cube of 
Y, from the decay of the explicit cluster propagators a large power of the product for all tubes of 
the factorial of the number of the vertices produced by the explicit cluster functional derivations 
and integration by parts [R] . This is useful because we have no symmetry factor in v ^y for these 
vertices. Joining this fact to (IV.2), we have a factor v( ^y where v(A) is now the total number 
of vertices in the polymer localized in A, no matter whether they are hooked to explicit cluster 
derivatives or not. 
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For any field , we call k{ij)) its scale k and Z($) the scale of the localization tube of the 
vertex to which it is hooked. We call A;(^>) the localization tube of the vertex to which it is hooked, 
not to be confused with the tube A($>) of scale fc($>) in which it lies. Remark that A;($) C A($). 

By a combinatoric factor K\ = M^/ 4 ^*^) - ^ we can sum over all scales of all legs of any 
vertex. In what follows e is a very small number. 

By a combinatoric factor n 2 = M - - 7-3 *^") per vertex of localization scale k(v) we can integrate 
over the position of that vertex in its localization cube. 

Then at each level k let us consider the piece Y k made of the tubes of Y of scale k. Each tube 
contains a certain number /(A) of fields (or anti fields) of scale k lying in this cube. 

We decompose the covariance C k (h,T) as C k (h,T) = J2a A'ex k XAC k (h,T)xA' and corre- 
spondingly we expand the fields. Let us introduce the scaled distance between A and A' as 

d(A,A') = inf (M fc |x - x| + M j \r - f |) (IV.3) 

By a combinatoric factor « 3 = (1 + d(A, A')) 5 we can choose, for each field or anti field, the tube 
(of the same scale) to which it contracts. 

Let /(A, A') be the number of fields (respectively antifields) localized in A that contract to a 
field localized in A'. By Fermionic rules of integration, /(A, A') = /(A', A). We obtain therefore 
after Gaussian integration for each scale a product of determinants 

n dct (A,Ao (iv.4) 

(A,A')er 2 

where det(A,A') is a /(A, A') by /(A', A) = /(A, A') determinant of propagators C- 7 ' . Applying 
the Hadamard's bound | detail < ?W 2 sup |ay|™ we obtain by (II. 5a) for these determinants a 
factor /(A, A') /(A ' A ' )/2 (1 + d(A, A'))-p/(a,a') m O+/c)/(a,a')^ Ta ki n g p larger than 10, the factor 
(1 + d(A, A'))~ p / 2 per field absorbs the combinatoric factor k 3 . Joining together the lines power 
counting JJ k Y\^ A A<)eY 2 M^ +k ^^ A ' A ) given by the determinant and by the explicit propagators to 

the integration factor k 2 we obtain the power counting factor k ^ 2 ^=1,2,3^*' k ^ per ver- 
tex, which is bounded by a factor M^'-^^ 1 /^) per vertex times a factor M ( ^/ 4 1+ e ^ k ^- 1 ^ 
per field. Therefore combining this factor to the factor JJ A v( ^y we have a factor 

M ((3/4) + e)(K^ , )-^ , ))„ (A ^ ) )) -l/4 (IV . 5 ) 

per field. After having absorbed the combinatoric factor Ki we also have still a factor M^ - ^")^ 1 / 4-2 ^ 
per vertex. 

Now we bound the factorials of Hadamard's inequality using in particular /(A, A') = /(A', A): 

n n /(a,a') /(a ^ a ' )/2 <n n /(A^o^^/^Ay^')/ 4 

k (A,A')eY 2 k (A,A')eY 2 

< n n (/(a, AO + /(a, a')) /(a ' a ' )/4+/(a ' a ' )/4 

fc (A,A')e^ 2 

<J] J] /(A)^A' /(A,A')/4+/(A,A')/4 

fc Aen 
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< n n /(a) /(a)/4 < n n ^ /(a) /(a)!! i/2 av.6) 

k AeYfc k AeY k 

So up to an inessential constant per field, Hadamard's bound gives a product over all tubes of 
the sqare root of the number of all Wick contractions between all fields localized in this tube, not 
taking into account whether they are fields or anti-fields. (Remark that /(A) can be odd, in which 
case we put by convention /(A)!! = (/(A) — 2). ..5. 3.1, hence wc do not contract one of the legs). 

Then we have: 
Lemma IV. 1 




(IV.7) 



Proof : We can contract the legs starting with the ones which have largest value of the factor 
M ( 3 + £ /4)(i(^)-*(V')) t ,(A I (^)) ) and with a factor 

M( 3+e / 4 )( | ^')-fc(^')) 4u (A J (^)) 



< Y M( 3 + £ / 4 )( ; (V')-fc(^))4w(A ; (i/;))M( 3 + e / 4 )( i (^ )~ k $ ))4u(Aj(^')) (!V.8) 

(-)' (-)' 

we can choose to which leg it contracts, by choosing first with M^ 3+£ / 4 ^'^ ) _fc (V )) the localization 
tube of the leg to which it contracts (this is because there are 

M 3(k-k') tubeg f sca i e k 

in a tube 

of scale k' < k), then by a factor v(A/ (?//)) the vertex in this localization tube, then with a factor 
4 the particular leg. Collecting these factors proves the lemma. ■ 



Comparing the factor of (IV.7) to the one of (IV. 5) we see that we get a constant per field. 

It remains to sum over the positions of the tubes in Y. This is standard ordinary power 
counting and it can be done with the decay of the cluster propagators in the horizontal direction 
(at fixed k) plus the power counting of the vertices involved in the vertical connections, using the 
fact that no connected subgraph has less than two external legs and less than three vertices. 

Finally the constant per field is controlled by the smallncss of A. This completes the proof of 
Theorem III. 1 . Taking into account the global fixed factor M J ( 4 ~ 2p ) due to our norm, Theorem I 
follows. 

Appendix I 

In this appendix we explain the slight modifications to adapt the proof of Theorem I to all 
space dimensions greater than three. 

The Fourier transform of the S function of the sphere S ,d_1 in TR d behaves as 1/ |a;| 4 )/ 2 at 
large x. Therefore reslicing the corresponding propagator in auxiliary slices C J ' fe we have 

|C J '' fe (£,0| < K{p)M j+k{d - 1 ^ 2 {l + M fe |x - x|)- p (l + M j \r - f\)- p (A.I.I) 

The power counting of a graph with n vertices at scale k and e = 2p external legs at scale 
j is M^ d+1 ^ d+1 ^ 2 ^ M^- k ^ n+ P^ d '^/ 2 - d K This is obtained by combining the power counting 

M (j+k(d-l)/2)(2n-p) of the legg at scale fcj and the powcr counting M -(j+dk){n-l) q{ t he Tl - 1 

vertices integration. Therefore the auxiliary problem is superrenormalizable in any dimension. 
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There is only a finite number of graphs to renormalize, those satisfying n + p(d — l)/2 < d, which 
are always vacuum (p = 0) or two point (p = 1) subgraphs. The main result of this paper, 
Theorem I extends therefore without difficulty to any dimension, with the prefactor M^ 2p ^ in 
(1.15) replaced by M( d+1 -P( d+1 )/ 2 W. 

ft is interesting to remark that as d — > oo more and more 2-point subgraphs diverge, and 
we obtain real Fcrmionic models (in integer dimensions) whose power counting mimick the ill- 
defined (p 4 models in fractional dimensions d — > 4. The "critical" dimension at which our problem 
becomes just renormalizable is d = oo, and it would presumably be interesting to set up for it a 
\jd expansion which would be the analogue of the Fisher- Wilson 1/e expansion for (j) 4 models. 

Appendix II 

In this appendix we include a discussion of why the sector decomposition used in [FMRT1] in 
the two-dimensional case does not work in d = 3. Since the volume of the j-th momentum shell is 
large, it is natural to decompose ths shell into M - ^ -1 ^ sectors of side AP, to get fields which arc 
supported on cells of unit volume in phase space. This canbe done through a smooth partition of 
unity 

M -(d-lW 

1= ^™(P)' Vm{p) = Vm f prfciT 

m=l VIPI 

of the Fermi surface, where r\ m is supported on the union of the m th sector, S m , and its neighbors, 
whose number is at most 3 d_1 — 1. There is a corresponding decomposition of the covariance 



M 



-(d-l)j 



m— 1 



where 



and of the fields 



J (2TT) d + L ip - e(p) 

M -(d-l)j M -(.d-l)j 

^ = i> u ' m \ ft = Y ft J,m) - 

m— 1 m— 1 

The standard power counting bound for an individual graph is still easy when there are 
sectors. First, one selects a spanning tree for the graph. To each line not in the tree there is 
a corresponding momentum loop, obtained by joining its ends through a path in the tree. This 
construction produces a complete set of independent loops. Ignoring unimportant constants, each 
propagator is bounded by its supremum M - - 7 . The volume of integration for each loop is now 
M( d+1 W. A priori, there is one sector sum with M^^ 1 ^ terms for each line. But, by conservation 
of momentum, there is only one sector sum per loop. Thus, if there are n vertices and e = 2p 
external lines, the supremum in momentum space of the graph is bounded by 

Yl M~ j Y[ M (d+1) - 7 'M _(d_1)j = M^ (4rl_2p)/2 M 2j ' [(4 ™ _2p)/2_(n_1)1 

lines loops 

= M5( 4 - 2 P)J. 



In the course of a non-perturbative construction, estimates cannot be made graph by graph 
because there are too many of them. Rather, the perturbation series must be blocked and the 
blocks estimated as units. The blocks are estimated using the exclusion principle to implement 
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strong cancellations between the roughly n! 2 graphs of order n. However, once the series is blocked, 
momentum loops can't be defined and the argument leading to the estimate above cannot be made. 
Conservation of momentum has to be implemented at each vertex, rather than through loops. Even 
though the volume cutoff A breaks exact conservation of momentum, many of the M -2 ^ -1 ^ terms 
in the sector sums for a general 2i?-legged vertex must be zero. 

In [FMRT1] the following lemma is proved: 
Lemma 2 Fix m e Z d+1 and I > 2 . Then, the number of 2l-tuples 

{Si, ■ ■ ■ S21} 

of sectors for which there exist ki G IR'', i = 1, • • • , 21 satisfying 

k^eSi, \ki -m < const M j , i = 1, • • • , 21 

and 

|ki + ••• +k 2t \ < const (1 + \m\)M 3 

is bounded by 

const ^(1 + |m|) d M" (2 ^" 1)(d " 1)j 'M j {1 + \j\8 d , 2 8 ta } . 
In particular, for a four legged vertex, the number of 4-tuples is at most 

const (1 + \n\) d M(- 3d +^ {1 + \j\8 da } . 

Here, k' = ^ denotes the projection of k onto the Fermi surface. 

This lemma although instructive, is neither for d = 2 nor for d > 2 powerful enough for 
the non-perturbative construction. However for d = 2 the number of active sector 4-tuples at a 
vertex is of order \j\M~ 2 i ; ordinary power counting at a vertex can accomodate at most M~ 2j , 
hence there is only a "logarithmic" power counting deficit of order \j\ at each vertex. For d = 3 
the number of active sector 4-tuples at a vertex is oforder M~ 5j , and ordinary power counting 
at a vertex can accomodate only M~ 4j so that the deficit at each vertex is a full power M~K 
This deficit is worse and worse as d increases beyond 3. In three dimension this difficulty can be 
geometrically thought as the twist or torsion which occurs in the figure made by four momenta of 
unit length adding up to (this figure is no longer a planar rhombus as in two dimensions). 

To circumvent the extra logarithm at d = 2, in [FMRT1] we divided the Fermi surface into 
sectors of length M^ 2 rather than M J . Taking into account the anisotropic spatial decay associated 
to these new "rectangular" sectors beats the logarithmic deficit. If we apply the same idea for d = 3, 
the number of active 4-tuples sectors at a vertex becomes of order M' b ^ 2 but the power counting 
in the manner of [FMRT1] gives M~ 2j , so that there still remains a deficit of order M~^ 2 per 
vertex. Making the sectors still longer does not work because the anisotropic decay no longer 
holds since the Fermi surface does not look flat for lengths bigger than M~H 2 . Remark that this 
difficulty is presumably linked to the graph G2, which as we recall is in M~ J A 2 in a cube of size 
M-i. 

The conclusion is that even the single slice theory in d > 3 is a non-trivial theory that contains 
a rather non-trivial renormalization, like <j)\. It is this renormalization that can be analyzed by 
means of the auxiliary scales used in this paper. 

Appendix III 

In this appendix we discuss briefly the present situation for the full (multi-scale) theory. 
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The most straightforward strategy in the case of multi-scale Fermionic models consists in 
writing the effect of inductive integration of higher scales in the form of an effective action for the 
single scale model, and to check some inductive estimate for this effective action. This is the road 
followed in [FMRT1]. However when combined with the Hadamard estimate on determinants, this 
approach runs into a major difficulty: even with a quartic bare action the effective action contains 
irrelevant operators of degree 6 and more in the fermionic fields. The Hadamard inequality when 
applied to such operators does not lead to a convergent answer (it develops divergent factorials). 

The obvious thing to do seems to treat the theory as a bosonic one. For bosonic theories 
we know well that one cannot simply exponentiate all the effective action, as in Wilson's original 
program, but that one must keep the irrelevant effects coming from higher scales in the form of a 
polymer gas with hardcore constraints [B] [R] . A typical stability estimate in this strategy is to prove 
that the "completely convergent" multi-scale polymers (i.e. those who do not contain any two or 
four point subcontributions) form a geometric convergent series. At the moment our best estimates 
with the Hadamard method indicate that the sum over such "completely convergent" polymer, 
extending over at most N different scales, converges provided the (bare) coupling constant A is 
small, and the number of scales is limited to N < e/ X (e being a small number). This is encouraging, 
but tantalizing. Indeed although this last estimate is much better than what a sector analysis would 
provide (namely N < K\ log A|), it is still not sufficient. For the construction of the BCS3 theory, 
where the BCS gap cuts the renormalization group flow at a scale Abcs — M~ N ~ e -0 ' 1 ^, 
we would need N ~ K/X (where K is not particularly small). For the construction of non- 
superconducting Fermi liquids as in [FKLT], the renormalization group flows for ever and we 
would need N independent of A. Therefore at the moment we consider that the question of the 
non-perturbative stability of these theories remains open, but we hope that the Hadamard method, 
clearly much better than the sector method, is a step towards the final solution. 
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